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We study the (exclusive) kr jet algorithm using effective field theory techniques. Regularizing 
the virtualities and rapidities of graphs in the soft-collinear effective theory (SCET), we are able 
to write the next-to-leading-order dijet cross section as the product of separate hard, jet, and soft 
contributions. We show how to reproduce the Sudakov form factor to next-to-leading logarithmic 
accuracy previously calculated by the coherent branching formalism. Our result only depends on 
the renormalization group evolution of the hard function, rather than on that of the hard and jet 
functions as is usual in SCET. We comment that regularizing rapidities is not necessary in this case. 



I. INTRODUCTION 

Jets are important for understanding the background 
to new physics being investigated at the Large Hadron 
Collider. Jet production is a multi-scale process that 
involves the large energy of the jet, Q, and its small in- 
variant mass, TOj 0t . A hierarchy of scales Q 3> mj e t gives 
rise to large logarithms of the form L = ln(Q 2 / m 2 ct ) 1 
in perturbative calculations. These logarithms manifest 
in the jet production rate in the form 



R — R nm a™L r 



(1) 



where a s is the strong coupling constant. Even when 
a s -C 1, the large logarithms will ruin perturbation the- 
ory when a s L 2 ~ 1. 

Well-known perturbative QCD (pQCD) techniques 
based on factorization theorems |lj and the coherent 
branching formalism [5] can sum these logarithms by 
writing the series ([!]) as 

R = C(a s )E(a s ,L), (2) 

where 

oo 

C{a s ) = C n a n s , (3) 

n=0 

ln£(a s ,L) = Lf (a s L)+ j x (a s L) + a s f 2 (a s L) + ... . 

The coefficient function C'(a s ) contains no large loga- 
rithms L, while £(cv s ,L) sums the logarithms. The /o 
term sums the leading logarithms (LL), the /i term 
sums the next-to- leading logarithms (NLL), and the f n >2 
terms sum the subleading logarithms. In this paper we 
will always refer to the logarithmic order in the exponent 
([3]) as opposed to the logarithmic order in the perturba- 
tive rate 0. 

An example of a jet definition is the (exclusive) k? jet 
algorithm [3J [3] , proposed to resolve the exponentiation 
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issue of the earlier JADE algorithm [31 [51 [5] • The kr and 
JADE algorithms combine final-state partons into jets 
using a distance measure y^ for all pairs of final-state 
partons {i, j}. If the smallest y^ is smaller than some 
pre-determined resolution parameter y c , then that pair of 
partons are combined and all the j/y's are re-calculated. 
The procedure is repeated until all j/y > y c , and these 
pseudo-partons are then called jets. The kx algorithm 
measure for e + e~ — > jets is 



yij = 2(1 - cos 6^) 



mm(E 2 ,E 2 



(4) 



where Q is the centre-of-mass energy, Oij the angle be- 
tween the final-state pair, and E^j their respective en- 
ergy. We are interested in a two-jet final state where the 
cut parameter is small. Jets in the y c <C 1 region have 
small mass mj Ct ~ \fycQ *C Q, which gives rise to large 
logarithms L = ln(l/y c ). The k^ dijet production rate 
has been calculated using the coherent branching formal- 
ism to full LL accuracy in 00] and partial NLL accuracy 
in [7] . Clustering effects among multiple gluon emissions 
generate unsummed logarithms that start at 0(a 2 L 2 ) in 
the exponent _S| and ruin the NLL summation of [7]. 

Effective field theory (EFT) techniques offer another 
approach to summing the large logarithms. Using EFTs 
has the advantage of using the renormalization group 
(RG) to sum the large logarithms, as well as providing 
a systematic approach to power corrections. In [9] the 
k-r dijet rate was calculated using soft-collinear effective 
theory (SCET) to next-to-leading order (NLO). SCET 
[TUhTS] describes QCD using highly boosted "collinear" 
fields and low energy "soft" fields. SCET has previously 
been successful in calculating jet shapes [Ml El, where it 
automatically separated the hard scattering interaction 
from the highly boosted interactions in the jets and from 
the soft radiation between them. Such a separation al- 
lows the rate to be written as the convolution of hard, 
jet (one for each of the dijets), and soft functions, 



R = H x J x J x S, 



(5) 



each of which depends on a different scale. These func- 
tions are then run individually to a common scale for 
logarithm summation. The authors of [9], however, were 
unable to use dimensional regularization to regulate the 
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individual NLO collinear and soft graphs of the kr dijet 
rate, making it unclear how to write the rate as separate 
jet and soft functions as in (J5J). 

Recently [TH1 US] a new regulator capable of regulating 
these divergences has been proposed. This new "rapidity 
regulator" effectively places a cut on the rapidities of the 
fields [in], enabling the rate to be written as separate 
scheme dependent jet and soft functions. The rapidity 
regulator was used to sum logarithms in the jet broaden- 
ing event shape [TH] [TS] [JJ5] , which has a similar issue at 
NLO to the fcr dijet rate. The introduction of the rapid- 
ity regulator opens up the possibility of the RG running 
in another scale v, in analogy to the usual RG running 
scale [i of dimensional regularization. 

We propose to extend the work of [9] using the new 
rapidity regulator and investigate how to write the kx 
dijet rate as the product of hard, jet, and soft functions 
as m (H). Our work provides another application of the 
rapidity regulator. As in [3] H] [7] we assume a factor- 
ization theorem, which allows us to interpret the SCET 
collinear and soft graphs as the jet and soft functions 
that are run using the RG. We can then use the RG to 
attempt to sum the large logarithms. We find that we re- 
produce the coherent branching formalism result [7] but 
that neither approach sums the logarithms generated by 
clustering effects [5]. A similar result was recently found 
for the inclusive fc^ algorithm [2D] . 

The summation of the logarithms in the dijet rate 
using SCET only requires the running of the hard func- 
tion to NLL accuracy. The jet and soft functions act as 
a single soft function S = J x J x S that reproduces the 
infrared physics of QCD and depends only on a single 
soft scale. For NLL accuracy, it is unnecessary to define 
separate scheme-dependent jet and soft functions using 
the rapidity regulator. 

The rest of the paper proceeds as follows: in Sec. [11] we 



review the NLO results and issues of [SJ, and in Sec. |III| 
we show how the rapidity regulator solves these issues. In 
Sec. [TV] we show our final result with NLL summation and 
compare with the coherent branching formalism result. 
We discuss the interpretation of our results and the utility 
of the rapidity regulator in Sec. [Vj We conclude in Sec. 



II. REVIEW OF PREVIOUS WORK 

The kr algorithm was previously studied using SCET 
in [9]. SCET is the appropriate EFT to describe QCD 
with highly boosted massless fields. Collinear fields de- 
scribe the boosted particles, and soft fields describe the 
low-energy particle exchanges. The interactions within 
each sector (soft, collinear in each direction) decouple 
from one another and are described by a copy of QCD 
[15] . The interactions between sectors in the full theory 
are reproduced in the currents via Wilson lines [TTHT5] . 

The appropriate SCET operator for dijet production 
where n and n are respectively the light-like directions of 
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(a) (b) 
FIG. 1. QCD diagrams for real emission 



the jets is [Hll2"T] 

2 = [Z n W n ] [Y^TY n ] 



(6) 



where £, n ,n is a two-component n- or n-collinear spinor. 
The Wilson lines are defined in momentum space as 
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(7) 



with Wn and Y n defined analogously. Here V 11 is the 
momentum operator that acts on the gluon fields. The 
fields A s ,A n , and A n represent soft, n-, and n-collinear 
gluon fields respectively. The matching between QCD 
and SCET is well known [22] and gives the matching 
coefficient 



C 2 (/x) = l + 



2ir 



1,2 C 2 
2 hl =G* 



and MS counterterm 
ci s Cf 



z 2 {n) = 1 + 



2tt 



3 



3 ^ 
2 ln =Q^ 



7T~ 

12 



(8) 
(9) 



The ellipses denote higher orders in a s . The matching 
coefficient reproduces the UV physics of QCD. 

The e + e _ — > 7* -> dijet rate is calculated in SCET by 
summing the collinear and soft diagrams and integrating 
over the appropriate phase space. Generally the rate is 
written in the form 



1 dcr dijet 



H x (J x J xS) 



(10) 



where CTo = (47ra 2 /Q 2 ) e'j is the Born cross sectior^j 
The soft contribution 5* describes the interaction of the 
soft fields, while the hard function H captures the physics 
of the hard initial interaction. The hard function is de- 
fined to be H = I C 2 1 2 - The jet contributions J and J 



For dijet rates via a Z°, only the Born cross section is modified. 
This is irrelevant for our calculation. 
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TABLE I. Phase space constraints for NLO real emission for 
the fcr algorithm. The constraints are plotted in Fig. [2] 



describe the interactions of the n- and n-collinear fields 
respectively. 

For perturbative calculations, the contributions in ( 10 ) 
are individually written as 



F(ji) = 1 + FW (/i) + {n) + . . . (11) 

where F = H,J,J,S and F^ is the 0«) term. The 
two QCD diagrams that contribute to real emission at 
NLO are shown in Fig. [l] In SCET, the gluon can either 
be soft, n-, or n-collinear, resulting in six graphs that 
must be summed. We write all momenta in lightcone 
coordinates = (n ■ q,n ■ q,qx) = (q + , q~ , q±). We 
adopt the convention of [5] and use the symbol k <C Q 
for soft momentum, and p ~ Q for large momentum. 
Contributions from the NLO collinear and soft graphs in 
dimensional regularization are given by integrating the 
corresponding differential cross sections over the relevant 
phase space PSf [S] 

PSs dktdK W^ 

dk^ dk^ 

'PS n 



5 (1) (M) = ^/. 



JW( M ) 



2tt 
2tt 



ft 





(1-e) 



Qk 3 

Q-Ps 



(13) 



ft 



2^r 



dk^ dk^ 
PSo (ktK) 1 



+ r 



(14) 



where J is referred to as the naive collinear graph and Jo 
the zero-bin. The "true" collinear contribution requires a 
zero-bin subtraction [23] and is denned as J(u) = J(/x) — 
Join). We have introduced / e = u 2£ e e7i5 /r(l - e) for 
later convenience. The rl-collinear graph is the same as 
the n-collinear graph at NLO, J (1) (") = J (1) (")- 

The relevant NLO phase space constraints in SCET 
are found by applying the kx measure Q to the qqg 
final state and expanding in k -C p, Q- At leading order 
in power counting the fermions must be collinear, and 
we define to be in the direction of the quark. The 
constraints for a soft and n-collinear gluon are shown in 
Table [T] and plotted in Fig. [2j The constraints for an n- 
collinear gluon are the same as those for an n-collinear 
gluon with "+" and "— " interchanged. 

In [5] it was found that the NLO soft graph can be 
written as 



7T er(l - e) \VcQ 



1 (1-^) 
dx { - ^ 



where the ellipses denote terms that are properly regu- 
lated in dimensional regularization. The integral in ( 15 ) 



is not regularized as x — > 0, and this means that inter- 
preting the soft function as the sum of the soft graphs as 
in ( |To| is not well defined. However, it was noted in [S] 
that the NLO zero-bin can be written as 

<CV ( u 2 V f 1 dx 



eT(l-c) \y c Q 



(16) 

where again the ellipses denote terms that are properly 
regularized. The x — > divergence in this integral is the 
same as the soft graph. Because Jq 1 ^ (/1) enters into both 
J(p) and J(/i) with a relative minus sign compared to the 
soft graph, the total rate (l/a )(da'^ ct /dy c ) is properly 
regularized at NLO as expected. 

Decomposing the rate as separately regularized jet and 
soft functions as in (10), where J and S are respectively 



(12) M 



the collinear and soft graphs, is therefore not possible 
using pure dimensional regularization. The issue of sepa- 
rately well-defined functions comes from how phase space 
is being divided between the collinear and soft graphs in 
this scheme. The soft graph is being integrated over the 
region k$ — > while keeping k^k^ < y c Q 2 . This is a 
highly boosted region, and is more naturally associated 
with the jet function than the soft function. The jet 
broadening rate has a similar issue in SCET as shown in 



As pointed out in [9], the soft graph can be regulated 
using a different scheme such as a cut-off regulator. The 
cut-off regulator removes the contribution of the afore- 
mentioned region from the soft graph and regulates the 
integral in ( 15 ). The jet broadening rate can also be reg- 



ularized using a cut-off. The cut-off regulator, however, 
is not very attractive as it is not gauge invariant, making 
it hard to run using the RG. It is also unclear how to 
define it in the naive collinear calculation. 

Another scheme also studied in [5] is to use offshell- 
ness as an infrared regulator, while using dimensional 
regularization to regulate the UV. Here, the small quark 
and anti-quark offshellness regulates the integrals in ( 15 ) 
and (16 1. However, the resulting collinear and soft con- 
tributions - including the virtual diagrams - are not indi- 
vidually infrared finite, even though these infrared diver- 
gences cancel in the total NLO rate as expected. There- 
fore it is again unclear how to interpret these as the jet 
and soft functions of (10 1. 



(15) 



In the next section we use the recently introduced ra- 
pidity regulator [T5J[T2] to separate the low energy theory 
into jet and soft functions associated with the collinear 
and soft fields respectively. 



III. NEXT-TO-LEADING-ORDER 
CALCULATION 

In this section we show how all the divergences in the 
phase space of the soft graphs are tamed with the intro- 
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FIG. 2. NLO &t dijet phase space for (a) n-collinear gluon (b) soft gluon and (c) zero-bin. These plots are taken from [9], and 
the bold arrows indicate that the plots extend to infinity. 



duction of the rapidity regulator [T51 QI5] . The rapidity 
regulator was used to solve the similar issue and sum 
the logarithms in jet broadening [HI [19]. The regulator 
acts as an energy cut-off in a similar way that dimen- 
sional regularization acts as a cut-off on the mass scale 
of loop momenta |24j . The form is similar to dimensional 
regularization and also maintains gauge invariance |19j . 
unlike a cut-off regulator. We will show in this section 
that using the rapidity regulator splits the NLO collincar 
and soft graphs into separately finite pieces. This allows 
us to interpret the jet and soft functions as the collincar 
and soft interactions respectively. 

The rapidity regulator modifies the momentum-space 
definition of the Wilson lines to [HI [19] 

o.\n-V\- 



w n =^2 ex p 
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exp 

perm 



v 



\ 2 p3\-r,/2 
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A n 



A, 



(17) 



Here V 3 pulls down the component of momentum in the 
spatial direction of the jet. The new parameter rj acts 
similarly to e in dimensional regularization. The param- 
eter w counts the number of emissions from a Wilson line, 
and is taken to one as rj — > 0. Implementing the rapidity 
regulator modifies the NLO collinear and soft graphs to 



2 w tf 

ol 8 Cf 



dk^ dk^ 



2n 



ps 3 fe + fc 3 ") 1+e 
(ktP3 ) 
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fe I dk£dp 3 



Qk 3 
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+2w 



2Q-P3 fPs 
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(18) 



a s C F 2 
2tt j 



dk^ dkr^ 
ps (ktK) 1+t V v 



Note that the phase space constraints PSf are not af- 
fected. The pure dimensional regularized functions are 
recovered in the 77—7-0 limit. 

Calculating the collinear and soft graphs is now 
straightforward. As has been previously demonstrated 
[19) . we must expand in r\ before e. As we are consider- 
ing the y c <C 1 region, all terms subleading in y c are also 
suppressed. 

The naive NLO collinear graph is 



a s C F 



2n 



Aw^ 1 



7T~ 

12 



In 2 



2+ ln2 
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w 2 (2 + lny c ) - - 



(19) 



We leave in w for now and will set it to one at the end. 
The logarithms cannot be minimized at any one scale be- 
cause we have not yet included the zero-bin subtraction. 
The NLO zero-bin contribution is 

12 9 „, r>2 



t(1)/ s a s C F 2 
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(20) 



Subtracting the zero-bin from the naive collinear graph 
gives the true (bare) collinear contribution 
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(21) 



The collinear logarithms can be minimized at fij = ^fy~ c Q 
and Vj = Q. The n-collinear contribution J(/i, v) is ex- 
actly the same as J(p, v) at this order in a s . 
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The NLO soft graph can be calculated similarly. The 
extra ry-dependent piece in ( 18 1 regulates the divergence 



of (15). The NLO (bare) soft graph is 



a,C 
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(22) 



In 



VcQ 2 



where the logarithms are minimized at the scales (is — 
\fy~cQ = v s- Note that the dimensional regularization 
scale of the soft graph is equal to that of the collinear 
graph, hj = us. 

Putting the collinear and soft graphs together, as well 
as the matching coefficient ^ and the counterterm 
the kj> dijet rate is 



1 daf jot 



O"0 



dy c 
a s C F 
2tt 



= H{(i)J{ii,v)J{n,v)S{n,v) 
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(23) 



In 2 y c 



3 In y c 
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which exactly reproduces the pQCD result El E] • All 
the graphs must be evaluated at the same (fi, v). Notice 
that the v dependence must cancel between the collinear 
and soft graphs because H is ^-independent. This is a 
general result and means that, when added together, the 
r\ dependence of the J, J and S counterterms must vanish 

eg. 

We find that unlike in [S], we can define the jet and 
soft functions in ( 10 1 as the collinear and soft interactions 
respectively. In the next section we show how to sum 
the logarithms using the RG by running each function 
individually. We then compare the summed expression 
to the coherent branching formalism result. 



IV. NEXT-TO-LEADING LOGARITHM 
SUMMATION 

We wish to calculate both / and /i of ^ to sum the 
logarithms and compare to [7] . Because we have two UV 
regulators, the jet and soft functions now run through a 
two-dimensional (fi, v) space. 

The renormalized function F is defined in terms of the 
bare function F B and counterterm Z F as F B = Z F F. 
Therefore, the anomalous dimensions in the two direc- 
tions of the (fi, v) space are found using 



7f(a*^) = 



_d_ 



In Zf 



(24) 



where F = H, J, S. The running of the coupling constant 
j3(a s ) = — 2a s e+0(a 2 ) is well-known and /3(w) — —r/w/2 



exactly [T5]. The counterterms of the jet and soft func- 
tions are found from (21) and (22) to be 
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(25) 



where we have set w — 1, /io 2 
here denote higher orders in a s 



= y c Q , and the ellipses 
The hard function coun- 



terterm Z H L = \Z 2 \ 2 = Z S Z 2 } as expected. The NLO 



anomalous dimensions are 
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(26) 



The hard anomalous dimension in the v direction van- 
ishes identically because Z 2 is ^-independent. For con- 
sistency in the running, we must have 



1h = 2jj +7s, 



(27) 



where 7f = (l F ,Y F ) = ~VlnZ F with V = 
From (26), we see that these conditions are satisfied at 
NLO. 

The anomalous dimensions allow the functions to be 
run to any scale. However, unlike in the usual case of 
only using dimensional regularization to regulate the UV, 
the hard, jet, and soft functions are now scalar functions 
defined over a two-dimensional (ji, v) space. Path in- 
dependence of running is equivalent to the curl of j F 
vanishing. This vanishing curl gives the condition 



M^-7f(m>2 / ) 
a/i 



^7f(m,^), 



(28) 



which, along with (27), must be satisfied to all orders in 
a s . We show in the Appendix that the soft v anomalous 
dimension can be written as 



7s (A*) = 7s (Mo) 



Mo M' 



^7sXM» 



(29) 



where the general form of the soft (i anomalous dimension 
is taken to be 

7s {^v) = T S [a a (//)] In ^ + 7s [a s (/«)]■ (30) 

Here Ts is called the cusp anomalous dimension. The 
7s(a*o) contains no logarithms and all the logarithmic 
dependency of 7^ (fi) is determined by the fi anomalous 
dimension. A similar expression to (29) appears in |19| . 
The hard anomalous dimension has a similar form as ( 30 ) 
with v = Q |16j . The jet anomalous dimension is com- 
pletely constrained by the hard and soft anomalous di- 



mensions from (27). 
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FIG. 3. Running each function from ((1f,Vf) to (fj,,f). 



We can use the above to solve the RG equations and 
sum the logarithms. Each function F(fi, v) must be 
evolved from the scale that minimizes its logarithms 
(ijlf, vf) to a common scale. The solution to the RG 
equations gives the running of each function 



(31) 



where we have chosen to run in v first but path indepen- 



dence is guaranteed with the use of (29). The summed 
rate for the path in Fig. [3] is therefore 



c dy c 



=H(/j, h )J(hj, vj)J(fij, fj)S(ns, vs) 
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(32) 

where we have run to a general (ju, v), and used the con- 
sistency equations (27) and path independence (29) to 



write everything in terms of the hard and soft running. 
Terms subleading to NLL accuracy have been suppressed. 
Because of path independence, we can choose any other 
path and get the same NLL terms. The summed rate is 
both fi- and ^-independent, as expected. 
The running kernels in ( 32 1 are defined as 
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where we denote r = a s (/xi)/a s (/i2). The coefficients T F 
and jp are given from the general form of the anomalous 
dimension ( 30 ) as 
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where from ( 26 1 we can read off 

= -8C F l° H = -12C F 

T% = 8C F 7s = 0. 



(35) 



The /3-function of the coupling constant a s also has an 
expansion 
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The two- loop running in the coupling constant a s (p) 
gives 
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The factor 
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(39) 



is the well-known ratio of the one- and two-loop cusp 
anomalous dimensions, K = Tp/T F [7J [TH], and is re- 
quired for the NLL summation. 

Choosing the scales that minimize the logarithms in 
the hard, jet, and soft functions 
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simplifies (32) to 



^- = ff(Q)J( Mo ,Q)J(Mo,Q)S(Mo,Mo) 

dy c 



(41) 



which sums the rate to NLL accuracy. From the above 
equation we can see that only the RG of the hard function 
is required for the summation to NLL accuracy. The 
action of running in rapidity cancels between the jet and 
soft functions. We will discuss this issue in more depth 
in the following section. 

We can now find the functions /o and f\ of ([3| from 



(41 ). The LL summation comes from setting 7^ = K = 



Px= 0. The NLL summation comes from the terms pro- 
portional to a single power of 7^, K, and Pi. Therefore, 
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Jo — -- 
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1 



m(l - x) 



(42) 
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FIG. 4. Plots of (411 for LL and NLL accuracy. NLO is the order in a s the coefficient function C(a s ) of |3| is taken to. The 
multipole expansion breaks down before L — > or y — } 1, and the inclusion of NLO terms in C(a s ) improves accuracy of the 
curve over this region. The procedure for calculating the error bands is described in the text. 



where a; = a s Cp/3oL/(4ir). Using (35) we see that the 



functions agree exactly with the coherent branching for- 
malism result [7]. We plot the summed rate in Fig. [I] as 
a function of ln(l/j/ c ). The maximum jet production is 
at around y c ~ 0.2, which corresponds to jets of mass 
y/y~ c M z ^ 40GeV for LEP - well above A QCD . 

The error in Fig.|4]is found by varying the scales Hh,j,s 
and vj^s in (32) by 2 and 1/2 of their values in (40). We 
vary the jet and soft scales together to maintain the /ij ss 
(is scaling. Varying the vp scales without varying fij 
produces no error due to the exponent ujs(^j — Mo, Mo) = 
0. We take a naive approach to estimate the correlated 
errors by varying fij and vp together, and taking the 
geometric mean of the resulting percent errors. 
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\ 71/ 
• • J 
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Q 



FIG. 5. Running the functions along a particular path. Note 
that, up to NLL accuracy, the summed result is independent 
of the path chosen and the final (/x, v) point. 



V. DISCUSSION 

That only the RG of the hard function is necessary 
for NLL accuracy suggests the kp dijet rate should be 
written as 



1 dat m 
o"o dy c 



(43) 



Here the new soft function 5(/i) = J(/i) J^S^) is the 
combined collinear and soft graphs and is well defined 
at NLO in pure dimensional regularization as seen in [pj 
and Sec. [TT] This new soft function is also infrared finite, 
as shown by using offshcllness to regulate the infrared di- 
vergences of the collinear and soft graphs [9j . By running 
the functions between fi^ = Q and /i^ = ^fy~ c Q = Ms. J 
the kp dijet rate (41) is reproduced to NLL accuracy. 



By choosing to run along the particular path in Fig. 
[5] it is clear that only the combined collinear and soft 
graphs are required for NLL summation. Along this path, 
the general form of the v anomalous dimension ( 29 ) be- 
comes 



Is (Mo) = ot s {n ) 



m>0 



IS a s 



(44) 



which contains no large logarithms. For N fc LL accuracy 
only the m < k terms are required. However, in general 
the 7g term, which is required for NLL accuracy, van- 
ishes as seen in the kp dijet rate above and all the cases 
in [15]. For N 2 LL accuracy, therefore, only the hard run- 
ning and the term are required. 
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FIG. 6. Running the soft function in a rectangular path back 
to itself. Naively this results in a LL phase when v 2 v\. 



VI. CONCLUSION 

We have studied the (exclusive) kr dijet rate using 
effective theory methods, and shown how to reproduce 
the coherent branching formalism result to NLL accu- 
racy. We must use the rapidity regulator if we wish to 
separate the NLO rate into regularized jet and soft func- 
tions. We have demonstrated how to sum to LL and 
NLL accuracy using the rapidity regulator in a path in- 
dependent way, which can be generalized to any process 
that has a factorization theorem. We comment that the 
rapidity regulator is unnecessary for summing the large 
logarithms to NLL accuracy in the example of the kj> 
dijet rate. The same accuracy can be achieved if we con- 
sider the combined jet and soft function and run to the 
common jet and soft scale. We also find that using SCET 
with a rapidity regulator does not account for clustering 
effects and cannot improve the coherent branching for- 
malism result. A more complicated SCET-like theory 
may be able to properly account for these clustering ef- 
fects, however, we do not explore such a theory in this 
paper. 
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APPENDIX 

Here we show, using the soft function as an example, 
how to obtain ( 29 ) , which allows us to sum to NLL accu- 



racy. Our argument relies on factorization of scales, the 



consistency condition (27), the vanishing curl (28 1, the 



general form of jh, and that the anomalous dimensions 
are defined perturbatively in a s . 

Factorization means that the anomalous dimensions of 
each function are sensitive only to scales relevant to it. 
Therefore, the v dependence of 7/ and 7s will only be of 
the form \n(y/Q) and ln(V//zo) respectively. The consis- 
tency condition ( |27[ ) requires that all v dependence of 7V 
and 7s must cancel to all orders in perturbation theory. 
As 7^ is cusp-like and 7^ vanishes, 7^ can have at most 
a linear dependence on /vp) and j F must have no v 
dependence. 

The appearance of ln(/x//xo) in 7V, on the other hand, 
is not constrained. These logarithms can show up in ar- 
bitrary powers, as long as they cancel one another in the 
sum 7s + 2-f j to reproduce jh- Fortunately these log- 
arithms have negligible effect on NLL summation. This 
fact is made clear by the particular path shown in Fig. 

t where these logarithms vanish in 7s (Mo)- Because of 
is and path independence, we suppress these terms in 

pi 



The fact that 7s is independent of v to all orders in a s 
is also fixed by the form of 7^ in ( 30 1 and the vanishing 
curl (28). Taking this general form and applying v(d/dv) 



to both sides of ( 28 ) yields 



(A.l) 



This means that v(d/di/)~/s(fj,, v) is independent of /x and 
in particular a s (n). Such terms do not exist in pertur- 
bation theory, unless 7s (a^) is independent of v. 

The full \i dependence of 7^ can therefore be obtained 
from 7^ via integrating (28): 



In (29) we choose fi' — /j,q such that all logarithms in 
7s O vanish and only the non- logarithmic terms remain. 
If (A. 2) is not used, then running the soft function in the 
closed path shown in Fig. [6] would result in a large phase 
that spoils the LL accuracy of the results when i/ 2 >■ v-y. 
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